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ABSTRACT 

We present a deep observation with the X-ray Multi-Mirror Mission of PSR B1823-13, a young pulsar with sim¬ 
ilar properties to the Vela pulsar. We detect two components to the X-ray emission associated with PSR B1823-13: 
an elongated core of extent 30" immediately surrounding the pulsar, embedded in a fainter, diffuse component of 
emission 5' in extent, seen only on the southern side of the pulsar. The pulsar itself is not detected, either as a point 
source or through its pulsations. Both components of the X-ray emission are well ht by a power law spectrum, 
with photon index E « 1.6 and X-ray luminosity (0.5-10 keV) Lx~9x 10^^ erg s“' for the core, and E Ri 2.3 and 
Lx ^3 X 10^^ erg s“' for the diffuse emission, for a distance of 4 kpc. We interpret both components of emission 
as corresponding to a pulsar wind nebula, which we designate G18.0-0.7. We argue that the core region represents 
the wind termination shock of this nebula, while the diffuse component indicates the shocked downstream wind. 
We propose that the asymmetric morphology of the diffuse emission with respect to the pulsar is the result of 
a reverse shock from an associated supernova remnant, which has compressed and distorted the pulsar-powered 
nebula. Such an interaction might be typical for pulsars at this stage in their evolution. The associated supernova 
remnant is not detected directly, most likely being too faint to be seen in existing X-ray and radio observations. 
Subject headings: ISM: individual (G18.0-0.7) — pulsars: individual (B1823-13) — stars: neutron — 
supernova remnants — X-rays: ISM 

1. INTRODUCTION 

Many of the energetic processes which take place in the af¬ 
termath of a supernova explosion are significant sources of X- 
ray emission. The supernova remnant (SNR) which results as 
the expanding debris collide with the interstellar medium (ISM) 
can generate both thermal X-rays from shock-heated gas and 
ejecta, as well as non-thermal X-ray synchrotron emission from 
the shock-acceleration of electrons. Meanwhile, the core col¬ 
lapse which triggers such explosions often forms a rapidly ro¬ 
tating central pulsar. Pulsars also produce X-rays, in the form 
of pulsed non-thermal emission from the pulsar magnetosphere, 
quasi-blackbody emission from the neutron star surface, or an 
extended pulsar wind nebula (PWN) powered by a relativistic 
particle outflow from the central source. 

Each of these processes dominates at different stages in the 
system’s evolution. Eor example, the youngest systems are 
generally well-suited for studying SNR blast waves and mag- 
netospheric emission from pulsars, while in middle-aged sys¬ 
tems one can better study the neutron star surface. Poten¬ 
tially the most challenging group of objects to understand are 
the “Vela-like” pulsars (Becker & Triimper 1997), which are 
at an evolutionary stage for which all the above processes 
may be taking place. Such pulsars, typified by the Vela pul¬ 
sar, have spin periods P ~ 100 ms, spin-down luminosities 
E = Att^IP/P^ ~ 10^®- 10^^ erg s“' and characteristic ages 
T = P/2P ~ 10-20 kyr (where P is the pulsar’s spin-period 
derivative and I is the star’s moment of inertia). Until recently, 
about 10 Vela-like pulsars were known; the highly successful 
Parkes Multibeam Pulsar Survey has approximately doubled 
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this number (Manchester et al. 2002). 

The Vela pulsar itself is sufficiently nearby (~ 250 pc) that 
its different emission processes can be clearly distinguished: X- 
rays from the neutron star surface, the pulsar magnetosphere, 
the surrounding PWN and the associated SNR have all now 
been identihed (Strickman, Harding, & de lager 1999; Lu & 
Aschenbach 2000; Helfand, Gotthelf, & Halpern 2001; Pavlov 
et al. 2001). However, the rest of the known Vela-like pulsars 
are at much greater distances, resulting in fainter, smaller and 
more heavily absorbed X-ray sources. It has thus been difficult 
to detect and interpret the emission from these objects. 

Radio observations of the pulsar B1823-13 have mea¬ 
sured a spin-period P = 101 ms and a period derivative P = 
7.5 X 10“'^ s s“', implying a spin-down luminosity E = 2.8 x 
10^® erg s“' and a characteristic age r = 21.4 kyr (Clifton et al. 
1992), all properties typical of a Vela-like pulsar. The distance 
to PSR B1823-13, as inferred from the dispersion of its ra¬ 
dio pulses, is 3.9 ±0.4 kpc (Cordes & Lazio 2002). Previous 
X-ray studies of PSR B1823-13 have been carried out with 
both ROSAT and ASCA. The ROSAT data showed a compli¬ 
cated source at the pulsar’s position, with a morphology best 
modeled as a point source embedded in a compact nebula of 
radius 20" and surrounded by an elongated diffuse region ~5' 
across (Einley, Srinivasan, & Park 1996). However, because 
of the low number of photons detected, it was not possible to 
determine whether the extended component of the source was 
a SNR or PWN, while the nature of emission from the pulsar 
itself also remained unclear. ASCA data also indicated the pres¬ 
ence of one or more compact components embedded in a more 
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extended region of emission, but lacked the spatial resolution 
needed to properly separate and interpret these sources (Saku- 
rai et al. 2001). 

As the first part of an effort to understand the X-ray emis¬ 
sion from Vela-like pulsars, here we present observations of 
pulsar B1823-13 with the X-ray Multi-Mirror Mission (XMM- 
Newton), the first observatory to provide sufficient angular, 
spectral and temporal resolution to separately identify all the 
processes described above. XMM-Newton also has much higher 
sensitivity (effective area 4650 cm^ at 1.5 keV) than previous 
missions, making it well-suited for studying faint sources such 
as seen here. 

2. OBSERVATIONS AND ANALYSIS 

Observations of PSR B1823-13 were carried out with XMM- 
Newton on 2001 October 16 and 2001 October 18, in two ob¬ 
servations each of length ~30 ks, as summarized in Table 1. 
The data described here correspond to the three X-ray imag¬ 
ing instruments on board XMM-Newton. The EPIC MOSl and 
MOS2 detectors were operated in the standard “full frame” 
mode, in which seven CCDs in each detector are used to pro¬ 
duce an approximately circular field-of-view of diameter 30', 
with a time resolution for each CCD frame of 2.6 s. The EPIC 
pn detector was operated in “small window” mode, in which 
only a central 4' x 4' square region is active, but for which the 
time resolution is 5.7 ms. In this mode a significant fraction of 
each 5.7-ms frame is used to read out the CCD, resulting in a 
dead-time fraction of 30% (Kuster et al. 1999). To avoid optical 
contamination of the field, the medium and thin blocking filters 
were used for the MOS and pn detectors, respectively. 

Initial processing of the data was carried out at the XMM- 
Newton Science Operations Centre (SOC). We analyzed the re¬ 
sulting event files using the XMM-Newton Software Analysis 
System (SAS), version 5.3.0. The data were first filtered to re¬ 
move hot pixels and other bad data, and to only allow standard 
event grades (patterns 0 to 12 for MOSl and MOS2, patterns 0 
to 4 for pn). Each of these data sets was examined for periods 
of high background by only considering events with energies 
in the range 10-15 keV. Times at which flares or high levels 
were seen in the X-ray count rate for this bandpass were ex¬ 
cluded. This latter filtering excluded a significant fraction of 
each observation; the resulting useful exposure times are listed 
in Table 1. 

2.1. Imaging 

For each observation and detector listed in Table 1, the data 
set was corrected was for vignetting losses using the SAS task 
EVIGWEIGHT. Next, an energy filter was applied to include 
only events falling in the energy range 0.5-10 keV, other ener¬ 
gies being dominated by background. Finally, all the MOS and 
pn data were combined to form a single image for each type 
of detector. For the MOS CCDs, the flat-fielded images have 
a plate scale of pixeE^; for the pn CCD, the flat-fielded 
image has a plate scale of 474 pixel"'. 

2.2. Spectroscopy 

The high background count-rate for XMM-Newton requires 
careful analysis in order to extract accurate source spectra. The 
background contribution to the data consists of two main com¬ 
ponents (see Lumb et al. 2002): a diffuse X-ray background, 

^http://xmm.vilspa.esa.es/ccf/epic/#background 
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assumed to have a smooth spatial distribution over the observed 
field, and a particle background, which shows spatial variations 
across the detector. 

For a source not much larger than the point spread function 
(PSF), we determine the spectrum by extracting data from both 
a small region enclosing the source, and from a reference spec¬ 
trum immediately adjacent. The spectrum of interest is then 
obtained by subtracting the reference spectrum from the source 
spectrum, scaling appropriately to account for the differing ar¬ 
eas of the extraction regions for the two fields. This approach 
assumes that because of the proximity of the two extraction re¬ 
gions, the differences in the background components are negli¬ 
gible between source and reference spectra. 

Analysis of extended sources is complicated by the spa¬ 
tially varying nature of the background. Both the diffuse X- 
ray background and the source photons are vignetted by the 
XMM-Newton mirrors, resulting in a background component 
that varies dramatically with detector position. The particle 
background, however, is unaffected by the mirrors, resulting in 
a background component that is largely independent of detector 
position. Given these complications, we follow the prescription 
given in Appendix A of Arnaud et al. (2002), in which refer¬ 
ence spectra from adjacent regions are used to correct for the 
X-ray background, and spectra from blank-field observations 
supplied by the SOC® are used to account for the particle back¬ 
ground. 

Once these corrections had been applied, spectra were re¬ 
grouped so that there were at least 50 counts per spectral 
bin for compact sources, and 100 counts per bin for ex¬ 
tended sources. For analysis of EPIC MOS data, we used 
the responses supplied by the SOC^ to provide information 
on the redistribution matrix and effective area of each CCD: 
specifically, for the MOSl CCD we used the response file 
ml_medv9ql9t5r5_all_15 . rsp, while for MOS2 we 
used m2_medv9ql 9t5r5_all_15 . rsp. For analysis of 
EPIC pn data, we generated our own response files using the 
SAS tasks RMFGEN and ARFGEN. Subsequent spectral fit¬ 
ting and analysis were carried out using XSPEC version 11.1.0. 

2.3. Timing 

In the “small window” mode used here, the EPIC pn CCD 
has sufficient time resolution to search for X-ray pulsations 
from PSR B1823-13. Due to improvements in the time-tagging 
of events since the time when pipeline processing was carried 
out on the data by the SOC, we completely reprocessed the pn 
data from the raw observation data files, using the SAS task 
EPCHAIN. After filtering the data as described in §2 above, we 
then corrected the arrival times of each event to correspond to 
a reference frame at the solar system barycenter. For a given 
extraction region and energy range, a pulsation analysis could 
then be carried out on the resulting arrival times. 

3. RESULTS 

3.1. Imaging 

The EPIC MOS image of PSR B1823-13 is shown in Fig¬ 
ures 1 and 2. Because a significant contribution to the back¬ 
ground emission is from high-energy particles, the vignetting 
correction which has been applied to the data gives the false 
impression that the background emission is stronger around the 
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edges of the fields. The O star HD 169727 (outside the extent of 
the region shown in Eigure 1) is detected as a prominent X-ray 
source; comparison of its position in our data with that in the 
Tycho Reference Catalogue (Hpg et al. 1998) suggests that the 
absolute astrometry of our observation is correct to within 0."5 
in each coordinate. 

The data clearly demonstrate that the X-ray emission associ¬ 
ated with PSR B1823-13 consists of at least two components. 
The brightest X-ray emission is seen immediately surrounding 
the position of the pulsar, in a compact region shown in Eig¬ 
ure 2 which we subsequently refer to as the “core”. However, 
the morphology of the core suggests that it is not an unresolved 
source, but rather is elongated along an extent of ^ 30" in 
the east-west direction. This extension can be seen indepen¬ 
dently but at lower resolution in the EPIC pn image (not shown 
here). We quantify this in Eigure 3, where we show the in¬ 
tensity profile of X-ray emission around the pulsar along both 
east-west and north-south axes. We compare these profiles to 
the XMM-Newton PSE for MOS1 at 1.5 keV, which we gen¬ 
erated using the King profile parameters included in the XMM- 
Newton calibration file XRT1_XPSF_0005 . CCF.* Above lev¬ 
els of ~ 25% of the peak flux in the image, the intensity profile 
in the north-south direction (i.e. the lower panel of Eig 3) ap¬ 
pears to be largely consistent with this PSE. However, the pro¬ 
file shown in the upper panel of Eigure 3 clearly demonstrates 
that along its east-west axis, the core region is significantly 
broader than the PSE, only falling to the apparent background 
level at offsets of « 15" on the eastern and western sides of the 
pulsar. Also present is a secondary clump of emission 20" east 
of the pulsar, seen at coordinates (J2000) RA 18*'26™14?5, Dec. 
-13°34'47" in Eigure 2, and and visible as a flux enhancement 
at an offset of -20" in the upper panel of Eigure 3. 

At larger angular scales, Eigure 1 demonstrates the exis¬ 
tence of a faint, amorphous component of emission extending 
to much larger distances from the pulsar than does the central 
core; we subsequently refer to this region as the “diffuse” com¬ 
ponent. This diffuse emission is located predominantly to the 
south of the core region. The morphology of the diffuse com¬ 
ponent is approximately elliptical, centered ^ 1.'5 south of the 
pulsar, with a major axis of ~ 5' aligned north-south, and a 
minor axis of ^ 4'. In all directions, the diffuse component de¬ 
creases in brightness with increasing distance from the pulsar. 
To the south, east and west of the pulsar, the diffuse component 
fades into the background over several arcmin; to the north of 
the pulsar the brightness profile decreases much more rapidly. 
We demonstrate this in Eigure 4, where we show the bright¬ 
ness profile of the X-ray emission along the north-south axis 
of the source. As a comparison, we have extracted the same 
profile from the blank-field observations supplied by the SOC, 
and have normalized the profile for the blank-field data to that 
of the source by comparing count rates in an identical source- 
free region. As shown in Eigure 4, a comparison of profiles 
between the source field and the blank field clearly shows that 
X-ray emission from the source extends at least 300" to the 
south of PSR B1823-13, but extends no more than ~ 40" to 
the north of the pulsar. Note that the ranges of offset shown 
in Eigure 4 correspond to the extent of the central CCD of the 
detector. At larger offsets the background dominates, and in 
these regions the differences between the background levels in 
the source field and in the blank field make it impossible to de¬ 


termine whether there is any significant signal above the noise. 
Eor this reason we cannot place useful constraints on the pres¬ 
ence of a faint limb-brightened shell of emission at larger radii, 
as might be produced by an associated SNR (see further discus¬ 
sion in §4.3 below). 


3.2. Spectroscopy 

We have extracted spectra for the core from both the MOS 
and pn data, using a circular extraction region, centered on the 
pulsar and with a radius of 18" (as shown in Eigure 2). Al¬ 
though only ~75% of the energy from an unresolved or slightly 
extended source is encircled at this radius, this extraction region 
minimizes contamination from the surrounding diffuse compo¬ 
nent. The core region is sufficiently compact that we can correct 
for background using the simple process of subtracting a refer¬ 
ence region from an adjacent source-free region. In this case, 
the reference region was located ~ 3' distant from the core, 
and was placed so as to avoid any contribution from the diffuse 
component of the emission. 

Six spectra resulted, one for each of MOSl, MOS2 and pn, 
for each of the two observations listed in Table 1. We then fit¬ 
ted each of power-law and blackbody models to these six spec¬ 
tra simultaneously, modified by photoelectric absorption — the 
best-fit results are summarized in Table 3. We find that the core 
region is well fit by a power law, with an absorbing column of 
Ah = (1.0 ± 0.2) X 10^^ cm“^ and a photon index E = 1 -6^ ^. At 
a distance to the system of 4d4 kpc, the unabsorbed 0.5-10 keV 
luminosity implied by this fit is Lx ~ 9dl x 10^^ erg s“', where 
we have corrected for the ~ 25% of the flux from this source 
which presumably falls outside our extraction region. A black- 
body model gives a much poorer fit to the spectrum. 

The diffuse component occupies a much larger area on the 
detectors than does the core region, and indeed fills the en¬ 
tire field-of-view of the EPIC pn image. Because there is no 
source-free region on the detector, it is impossible to accu¬ 
rately account for the background of the diffuse emission for 
the EPIC pn. Instead, we only considered MOS data in this 
analysis. We extracted a spectrum from a circular region cen¬ 
tered on (J2000) RA 18''26”14?4, Dec. -13°36'14"8, with ra¬ 
dius 200", as shown in Eigure 1. We excluded from this region 
a 30"-radius region centered on the pulsar, and 20"-radius re¬ 
gions around the point sources seen at (J2000) RA 18*'26™03?5, 
Dec. -13°36'41"7 and RA 18''26”10?6, Dec. -13°39'08"4. 
As a reference spectrum, we considered all emission outside 
the main extraction region, but falling within a square of side 
9(5 centered on the pulsar. Both source and reference spec¬ 
tra were corrected for the particle background by subtraction of 
blank field data, as explained in §2.2. All four MOS spectra 
were then fit simultaneously to a variety of absorbed models 
(see Table 3 for details); a power law, thermal bremsstrahlung, 
a synchrotron spectrum with an exponential cut-off (model 
“SRCUT”), and three thermal plasma models (“RAYMOND”, 
“MEKAL” and “NEI”). In all six cases, we find a good match to 
the data with an absorbing column of Nh ^ 1.2 X 10^^ cm Eor 
the power law fit shown in Eigure 5, the photon index is E « 2.3, 
while the thermal models all imply temperatures kT ~ 3.5 keV. 
The cut-off synchrotron spectrum can be extrapolated to radio 
wavelengths, where the corresponding 1 GHz flux density im¬ 
plied for the source is 15 ± 1 mJy.^ 


^See http://xmm.vilspa.esa.es/docs/documents/CAL-SRN-0100-0-0.ps.gz for more information. 
‘>1 mJy = 10-2® Wm-2 jj^-l 
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3.3. Limits on a Point Source 

While it is clear that the core component of X-ray emission 
surrounding PSR B1823-13 is extended, an important consid¬ 
eration is how much of the X-ray emission detected could be 
from an unresolved source potentially corresponding to the pul¬ 
sar itself. We here consider the upper limits we can place on the 
temporal, morphological and spectral contributions of a central 
point source. 

3.3.1. Timing Analysis 

In principle, PSR B1823-13 may be directly detectable in 
these data through its pulsations. Because any emission from 
the pulsar is embedded in the surrounding extended compo¬ 
nent, pulsations will inevitably be contaminated by this ex¬ 
tended emission, and the signal-to-noise ratio of the pulses will 
thus depend on the extraction radius and energy range in which 
data are considered. We therefore searched for pulsations by 
considering events in five possible circular extraction regions, 
of radii 2"5, 5", 10", 15" and 20", each centered on the pulsar 
position. For each extraction radius, we considered photons in 
the energy ranges 0.5-2, 2-5, 5-10 and 0.5-10 keV. We then 
determined the phase of each photon by folding the data at the 
ephemeris listed in Table 2, which was kindly supplied to us 
by M. Kramer from pulsar timing observations at the Jodrell 
Bank Observatory. Because our two observations were only 
separated by two days, we combined data from both runs when 
searching for pulsations. 

With the phases corresponding to the known pulsar period in 
hand, we applied to the data the test (Buccheri et al. 1983), 
a binning-independent method of looking for the nth harmonic 
of a pulsating signal in sparsely-sampled data. In the absence 
of pulsations, the quantity will be distributed like with 2n 
degrees of freedom. For each of the extraction radii and energy 
ranges listed above, we have calculated Z^ for our observations 
for harmonics in the range 1 < n < 5. Amongst these data, the 
most significant statistic is Z| = 9.5, which occurs for an extrac¬ 
tion radius of 15" and an energy range 0.5-10 keV. The prob¬ 
ability of such a signal emerging by chance is 4.9%, which is 
consistent with chance occurrence given the number of differ¬ 
ent extraction radii and energy ranges searched. We therefore 
conclude that there is no significant pulsed signal in the data. 

We can estimate an upper limit on the amplitude of any pul¬ 
sations as follows. For a sinusoidal signal of pulsed fraction 
/ embedded in an event list containing N photons, we expect 
Zj = 0.5f^N (Feahy, Eisner, & Weisskopf 1983). Since we 
do not know the optimal extraction radius in which pulses are 
most likely to dominate over surrounding extended emission, 
we conservatively consider a 20" extraction radius. This cor¬ 
responds to ~ 75% of the energy of a central point source'®, 
but also most likely suffers from significant contamination from 
surrounding emission. The resulting values of Zj can be used 
to determine upper limits on / for the entire core region. For 
0.5-2, 2-5, 5-10 and 0.5-10 keV, we determine upper limits on 
/ for a sinusoidal pulsed signal from PSR B1823-13 of 11%, 
4%, 14% and 2% respectively. The upper limits on the pulsed 
fraction for an unresolved pulse are a factor of two lower in 
each case. 

3.3.2. Imaging Analysis 

The intensity profiles in Figure 3 show that while cuts 
through the pulsar along a north-south axis are consistent with 

'®See http://heasarc.gsfc.nasa.gov/docs/xmm/uhb/XMM_UHB/nodel7.html 


a point source, those in an east-west direction demonstrate the 
presence of extended emission immediately surrounding the 
pulsar. However, unlike other cases in which an unresolved 
source is clearly superposed on a smooth nebula (e.g. Murray 
et al. 2002), there is no obvious decomposition of the profile 
seen in the upper panel of Figure 3 into two such separate com¬ 
ponents. The profile shown could be interpreted as a centrally 
peaked nebula, with no central point source, but it could equally 
be regarded as an unresolved source embedded in lower level 
nebular emission. Thus we must adopt a conservative upper 
limit on a point source, in which we assume that all the flux 
in the central bin corresponds to unresolved emission. Inte¬ 
grating under the appropriately normalized PSF, we find that in 
the energy range 0.5-10 keV, ^ 1010 counts in the combined 
MOS1-I-MOS2 data set could be from an unresolved source at 
the pulsar position. (Note that only ~ 75% of these counts 
would fall within the 18" extraction radius shown in Figure 2, 
owing to the large wings of the XMM-Newton PSF.) The cor¬ 
responding upper limit on the EPIC MOS count-rate from the 
pulsar is 0.0123 cts s“*. 

3.3.3. Spectral Analysis 

The spectrum of the core region is well fit by a power law. If 
the emission from the pulsar has a blackbody spectrum, then it 
should be spectrally distinct from the rest of the core. Thus, we 
can constrain the temperature of a central source by determin¬ 
ing what contribution a thermal source could make to the core’s 
spectrum (Slane, Helfand, & Murray 2002). Specifically, we 
adopt a hard upper limit on the foreground column density of 
Nh < 2 X 10^^ cm“^, above all values of Nh listed in Table 3. 
We assume that the neutron star’s radius as viewed by a distant 
observer is Roo =12 km, and that the distance to the source is 
4 kpc. We can then fit a two component spectral model to the 
core’s spectrum (as measured by the combination of MOSl, 
MOS2 and pn), in which we allow the photon index and flux 
of the power-law component to be free parameters, and find the 
maximum temperature of the blackbody component which still 
yields an acceptable fit to the spectrum. This results in an upper 
limit on the surface temperature of the blackbody (as viewed at 
infinity) < 147 eV, and an unabsorbed bolometric luminos¬ 
ity < 8.7 X 10^® erg s“'. We note that a blackbody model 
will overestimate a neutron star’s temperature when compared 
to more realistic atmosphere models (Pavlov & Zavlin 2000), 
strengthening the upper limits made from blackbody fits here. 

4. DISCUSSION 

4.1. Comparison with Earlier Results 

ROSAT observations of this source suggested the presence of 
three components of emission; an unresolved source centered 
on the pulsar, embedded in a compact region of radius 20", and 
further surrounded by a diffuse nebula of extent ^ 5' (Finley, 
Srinivasan, & Park 1996). From the ROSAT data, there was the 
suggestion that this latter diffuse component had a cometary 
morphology, with a fan-like tail extending out behind the pul¬ 
sar at a position angle ~ 210° (north through east). In the data 
presented here, we confirm the presence of both a compact re¬ 
gion surrounding the pulsar, and of a diffuse nebula at larger 
extents. However, our results differ from those of Finley et al. 
(1996) in two important respects. 

First, based on the analyses in §3.3 above, we can only put 
upper limits on the contribution of a central unresolved source. 



GAENSLER ET AL. 


5 


Given the many more counts detected here, we ascribe the ap¬ 
parent detection of such a source by ROSAT to statistical fluctu¬ 
ations in the presence of low signal-to-noise. We can estimate 
the expected contribution from the pulsar as follows. The X-ray 
spectrum of the Vela pulsar can be modeled as a blackbody with 
Roo = 2.1 km, Too = 128 eV and = 1.5 x 10^^ erg s“', com¬ 
bined with a power law component of photon index E = 2.7 and 
an unabsorbed 0.2-8 keV luminosity Lx = 4.2 X 10^' erg s ^ 
(Pavlov et al. 2001). The pulsed fraction in the 0.1-10 keV 
band is 7% (Helfand, Gotthelf, &. Halpern 2001). At a dis¬ 
tance of 4 kpc, and adopting Nh = 1 x 10^^ cm“^ as repre¬ 
sentative of the minimum of the range of possible absorbing 
columns derived in Table 3, we would expect to detect ~ 35 
EPIC pn counts in our observation from such a source. When 
embedded in the surrounding core region, the effective pulsed 
fraction would be < 1%, well below the limits determined in 
§3.3.1. The resulting MOS count-rate for such a source would 
be ~ 4 X 10 cts s ', also well below the upper limit on the 
point-source count-rate determined in §3.3.2. Einally, the ex¬ 
pected temperature and luminosity of the blackbody component 
are both much less than the spectral limits on such a source de¬ 
termined in §3.3.3. Thus the failure to detect PSR B1823-13 
in our data is entirely consistent with the emission properties 
expected of a Vela-like pulsar. Higher resolution observations 
with Chandra will be needed to separate out any unresolved 
source from the surrounding extended component, and to thus 
directly detect the pulsar. 

Second, while we confirm the presence of diffuse emis¬ 
sion surrounding the core region, we find no evidence for the 
cometary morphology proposed by Einleyetal. (1996). Rather, 
it seems likely this appearance resulted from the much lower 
sensitivity of their data (~ 280 total source counts for ROSAT 
PSPC compared to ^ 4200 for XMM-Newton EPIC MOS), 
and simply corresponded to the brightest regions of the larger 
source seen here. 

Our measurements also differ somewhat from the ASCA re¬ 
sults presented by Sakurai et al. (2001). While Sakurai et al. 
(2001) concluded that the X-ray emission consists of a core re¬ 
gion embedded in a more diffuse nebula, the diffuse compo¬ 
nent in their data seems symmetrically distributed around the 
pulsar, as opposed to the distinctly one-sided morphology seen 
here. Sakurai et al. (2001) find that the core and diffuse com¬ 
ponents can both be fit with power-law spectra, with photon 
indices of r= 1.9 ±0.2 and F = 2.0 ±0.1 respectively. While 
these spectra are broadly consistent with our measurements, the 
unabsorbed flux densities (2-10 keV) for the core and diffuse 
regions inferred from ASCA are ~ 2 x 10“'^ erg cm“^ s“' and 

2 X 10“'^ erg cm“^ s“' respectively, an order of magnitude 
higher than those inferred here. Since Sakurai et al. (2001) 
provide no information on their analysis procedure, fitted ab¬ 
sorbing column or assumed distance, it is difficult to assess 
whether there is a genuine discrepancy between our results and 
theirs. Pending a more detailed analysis of the ASCA data, we 
ascribe these differences to the much poorer angular resolution 
and sensitivity of ASCA when compared to XMM-Newton. 

4.2. Morphology of the X-ray Emission 

Spectral fits to the the diffuse component of the X-ray emis¬ 
sion cannot distinguish between thermal and non-thermal mod¬ 
els. In the case of a thermal spectrum, this emission could be 
interpreted as shock-heated gas from a SNR associated with the 
pulsar. However, the temperatures inferred by the correspond¬ 
ing fits listed in Table 3 are higher than seen even in very young 


(< 2000 yr) SNRs (e.g. Blanton & Helfand 1996; Gotthelf & 
Vasisht 1997). Pits with enhanced abundances can lower the 
inferred temperature (e.g Pivovaroff et al. 2001) — we have at¬ 
tempted to model the data in this way here, but the quality of 
the fit then becomes unacceptable. We thus conclude that the 
spectrum of the diffuse component is most likely non-thermal, 
and is best described by either the power-law or synchrotron 
cut-off spectral model. 

The synchrotron cut-off spectrum (model “SRCUT” in 
Table 3) corresponds to synchrotron emission from shock- 
accelerated electrons in supernova remnants (Reynolds & Keo- 
hane 1999). It is not unreasonable to suppose that a SNR asso¬ 
ciated with PSR B1823-13 is accelerating such electrons, and 
generating the diffuse component of the emission. Indeed, the 
photon index of a power law fit to the spectrum, F « 2.3, is sim¬ 
ilar to that seen from other SNRs emitting synchrotron X-rays 
when fit with such a model (Slane et al. 1999, 2001). Fur¬ 
thermore, the inferred 1 GHz flux density from the SRCUT 
fit, 15 ± 1 mJy, is well below the upper limit on 1.4-GHz ra¬ 
dio emission from this source (~ 50 mJy, 3a) as determined 
by Braun et al. (1989) from their radio observations of this 
field. However, in this interpretation the morphology of the X- 
ray emission should take the form of an approximately circu¬ 
lar limb-brightened region associated with the expanding SNR 
shock (e.g. Koyama et al. 1995). The amorphous, centrally- 
brightened emission seen here does not meet this expectation. 

We therefore interpret the diffuse component of the X-rays 
seen surrounding PSR B1823-13 as synchrotron emission pro¬ 
duced by a pulsar wind nebula, which we designate G18.0-0.7 
based on its Galactic coordinates. The fact that the diffuse com¬ 
ponent steadily fades with increasing distance from the pulsar 
supports the interpretation that PSR B1823-13 is the central 
energy source for this emission. The inferred photon index for 
the diffuse component, F ~ 2.3, is typical of PWNe (Becker 
& Triimper 1997; Possenti et al. 2002). The unabsorbed 0.5- 
10 keV luminosity for this source is Lx ~ 3dl X 10^^ erg s 
the corresponding efficiency with which the pulsar’s spin-down 
is converted into nebular X-rays is Lx/E = 1-1 X 10 a value 
also typical for young pulsars (Becker & Triimper 1997). 

Contained completely within the diffuse component, the core 
region is a compact region with a hard power-law spectrum cen¬ 
tered on the pulsar. Bright central structures are seen around 
many other young pulsars which power PWNe (e.g. Helfand, 
Gotthelf, & Halpern 2001; Gaensler et al. 2002a), and simi¬ 
larly we argue here that the core region is a central region of 
synchrotron emission powered by the pulsar. We thus inter¬ 
pret the core and diffuse regions of X-ray emission as both 
parts of a single PWN, of total unabsorbed luminosity (0.5- 
10keV)Lx«4t// X 10^^ erg s ^ and foreground column density 
Ah ~ 1.2 X 10^^ cm“^. Spectral fits with the “SRCUT” model 
demonstrate that a continuous synchrotron spectrum for G18.0- 
0.7 extrapolated to radio frequencies results in a 1-GHz flux 
density well below the sensitivity of radio searches for a nebula 
of this size, accounting for the non-detection of any PWN sur¬ 
rounding PSR B1823-13 at radio wavelengths (Braun, Goss, & 
Lyne 1989; Gaensler et al. 2000). 

In such a source, we expect that the pulsar wind flows freely 
out to the point at which the wind pressure is balanced by that 
from the nebula. In this region the wind decelerates, a ter¬ 
mination shock is formed, and synchrotron emission is gen¬ 
erated. Recent high-resolution observations with the Chandra 
X-ray Observatory have provided beautiful images of termina- 
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tion shocks around young pulsars. In some cases these images 
demonstrate that the termination region is not a spherical shock, 
but rather has a toroidal morphology presumably corresponding 
to a wind focused into the pulsar’s equatorial plane (Gaensler 
et al. 2002a; Lu et al. 2002). 

While we generally expect any X-ray PWN to be brightest 
near the pulsar, Figure 3 demonstrates that the core is a distinct, 
well-defined region superimposed on the much fainter diffuse 
background. We thus argue that the core component seen here 
represents the pulsar termination shock, its elongated morphol¬ 
ogy possibly indicative of anisotropic outflow from the pulsar. 
The radius of the termination shock is then « 15" = 03dA pc, 
an extent similar to that seen around the Crab pulsar and 
PSR B1509-58 (Weisskopf et al. 2000; Gaensler et al. 2002a), 
but about an order of magnitude larger than that inferred for the 
Vela pulsar (Helfand, Gotthelf, & Halpern 2001). We expect 
particles emitting in this region to have undergone minimal syn¬ 
chrotron losses, and to thus have a harder spectrum than the rest 
of the nebula (e.g. Slane et al. 2000; Gaensler et al. 2002a). We 
can directly compare the spectra of the core and diffuse compo¬ 
nents by fitting each to a power law with fixed Nh- The results 
of such fits are given in Table 3 for Nh = 1.0 x 10^^ cm“^ (cor¬ 
responding to the best-fit absorbing column for the core) and 
Nh = 1.4 X 10^^ cm“^ (corresponding to the best fit for the dif¬ 
fuse region). From these fits, we see that the core is indeed 
harder than the rest of the nebula at a statistically significant 
level, with a difference in photon index AF « 0.3-0.5. 

The X-ray luminosity of the core region is about 25% of that 
from the total nebula. This is somewhat higher than the ~ 10% 
contribution made by the central termination shock region in 
other such sources (Slane et al. 2000; Gaensler et al. 2002a; 
Fu et al. 2002), potentially indicating a difference in the effi¬ 
ciency or length scale over which electrons and positrons are 
accelerated in the shock seen here compared to other sources. 
However, we caution that a significant contribution to the flux 
in the core region could be from the pulsar itself (see §3.3.2 
above). 

If the core indeed represents the termination shock, the extent 
of this region then provides a constraint on the magnetic field 
in G18.0-0.7. Specifically, we expect pressure balance between 
the pulsar wind and the surrounding nebula. 


Pneb ‘ 


47rHr?c ’ 


( 1 ) 


where Pneb is the pressure in the PWN and H is the fraction of a 
sphere through which the wind flows. The limited angular res¬ 
olution of XMM-Newton makes it difficult to derive a specific 
geometry for outflow from the core morphology; we here leave 
H as a free parameter with a hard upper limit H < 1. The cor¬ 
responding nebular pressure needed to confine the pulsar wind 
is Pneb = 9.2Vr'^d1^ X 10“^^ erg cm“^. At distances far down¬ 
stream of the shock, we expect approximate equipartition be¬ 
tween particles and magnetic fields (Kennel & Coroniti 1984), 
so that Z?^/87r « Q.SPneb, where B is the mean nebular magnetic 
field. In this case we can infer B w pG. For elec¬ 

tron/positron pairs emitting at 1.5 keV (the peak of the X-ray 
flux seen in Figure 5), the corresponding synchrotron lifetime 

is tsynch ~ yr. 

The simplest expectation for a PWN is that, beyond the ter¬ 
mination shock, the emitting particles should flow out radially 
in all directions until they either radiate all their energy or reach 
the outer edge of the PWN. However, G18.0-0.7 differs from 


this in two important respects. Most fundamentally, rather than 
forming an approximate ellipsoid with the pulsar and termina¬ 
tion shock at its center (e.g. Slane et al. 2000, 2002), G18.0- 
0.7 is distinctly one-sided, showing little emission to the north 
of the pulsar but a broad diffuse region to the south. Also prob¬ 
lematic is the fact the southern region of the PWN extends to a 
distance L^5' = Ma, pc from the pulsar. If this extent represents 
the outer edge of the expanding nebula in this direction, then the 
expansion speed of the nebula is V w L/t w 3QQd4 km s“' (e.g. 
van der Swaluw et al. 2001). Alternatively, if the true nebular 
extent is much larger than that seen in X-rays, then we expect a 
flow speed at the edge of the X-ray nebula V w (rs/L)^(c/3) w 
250 km s“^ (Kennel & Coroniti 1984). Both these calculations 
are at odds with the speed needed to transport particles to the 
southern end of the PWN within their synchrotron lifetimes, 
L/tsynch ~ 600011“^/'*km s“^. Thus the size of the PWN 
is much larger than be easily accounted for if the nebula has 
expanded unimpeded on its southern side. 

One possible explanation for the PWN’s properties is if if 
the pulsar’s space velocity is supersonic with respect to the 
sound speed of its environment. In this case the PWN is con¬ 
fined through ram pressure and a bow-shock PWN results (e.g. 
Gaensler, Jones, & Stappers 2002). In such PWNe, the pul¬ 
sar is at one end of an elongated X-ray nebula; the PWN has 
a bright head at or near the pulsar position, with a progres¬ 
sively fainter tail extending opposite the direction of motion 
(e.g. Kaspi et al. 2001a). While there is no proper motion in¬ 
formation for PSR B1823-13 from which we might estimate its 
space velocity, the gross morphology of G18.0-0.7 could be in¬ 
terpreted as a tail extending to one side of the pulsar, suggestive 
of rapid motion of the central source in a northerly direction. 
However, there are at least two difficulties with this interpreta¬ 
tion. 

First, in other bow-shock PWNe, the shocked pulsar wind 
interior to the bow shock is tightly confined by the flow of 
surrounding material, and the X-ray tail correspondingly con¬ 
sists of a very narrow column of emission, marginally resolved 
or unresolved even with Chandra (Bucciantini 2002; Stappers 
et al. 2003). On the other hand, the diffuse emission seen here 
forms a very broad region behind the pulsar, showing no evi¬ 
dence for containment or collimation. Second, the length of the 
proposed tail region seen here, L « Ma, pc, is too long to be a 
trail of emission left behind by a fast pulsar. Even if the pulsar 
had an extreme space velocity VpsR = 2000 km s“', the lifetime 
computed above for X-ray emitting pairs limits the length of 

such a synchrotron “wake” to VpsR tsynch ~ pc, no¬ 

ticeably smaller than the measured extent of G18.0-0.7. Alter¬ 
natively, if the emitting particles have not been left behind by 
the pulsar, but are part of the wind flow driven by the pulsar, 
one is left with the uncomfortably high-flow speed discussed 
above, L/tsynch ~ 6000 km s“^. Wang et al. (1993) argue that 
such high-speed backflows behind the pulsar can be attained 
when the wind is forced to travel through a narrow nozzle, but 
no such collimated region behind the pulsar is seen here. 

A second possibility to explain the morphology and extent of 
this PWN is that the nebula is intrinsically symmetric around 
the pulsar, but that there is a relativistic bulk flow of particles 
within the nebula which results in Doppler boosting of the ap¬ 
proaching component over receding material. Such a situation 
has been proposed to explain the one-sided PWN seen around 
PSR B1509-58 (Gaensler et al. 2002a), in which a flow veloc¬ 
ity (3 = v/c 0.2 has been inferred. However, in PWNe for 
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which such high speed flows have been claimed, this outflow 
manifests itself in collimated jets of material directed along the 
pulsar spin-axis (Gaensler et al. 2002a; Lu et al. 2002) — the 
lack of any such collimated structures here argues against this 
possibility. The large extent of the PWN can be explained in 
this case, since the flow time from the pulsar to the edge of the 
nebula is only tfiow ~ 20/3“^ yrs (ignoring projection effects). 
However, at the relativistic speeds /3 ^ 0.1 needed for Doppler 
boosting to be effective, we And that tfio„ <C tsynch- We thus 
would expect minimal synchrotron cooling of the outer nebula 
(i.e. a photon index E « 1.6 as is seen in the core), in contradic¬ 
tion with the steeper photon index (F « 2.3) seen for the diffuse 
component of G18.0-0.7. 

A final alternative is that the system is in an evolved state, 
in which the reverse shock from a surrounding SNR has prop¬ 
agated back inwards, colliding with the PWN (Reynolds & 
Chevalier 1984). If different components of this reverse shock 
collide with the PWN at different times, emitting particles can 
be rapidly convected away from the pulsar, generating a neb¬ 
ula which can be significantly offset to one side of the pulsar 
(Chevalier 1998). Such an effect has been specifically pro¬ 
posed to account for the one-sided morphology of the PWN 
powered by the Vela pulsar (Blondin, Chevalier, & Erierson 
2001), and we propose that this process may similarly be oc¬ 
curring here for the Vela-like pulsar B1823-13. Such an asym¬ 
metric PWN will result only if the PWN crushing time-scale 
is of order or less than the lifetime of emitting pairs, and if 
the velocity of convective motions is sufficient to produce the 
observed PWN extent. Simulations of this process suggest 
a crushing time-scale ~ 2000 yrs and convective motions at 
speeds ~ 1000 km s“^(van der Swaluw et al. 2001). Assuming 
that the PWN had a radius of a few pc before interacting with 
the reverse shock, these values are roughly consistent with the 
observed nebular extent and the lifetimes of emitting electrons 
inferred above, and we thus think an asymmetric reverse-shock 
interaction the most likely explanation for the observed nebular 
morphology. One important consideration in this interpretation 
is that we do not directly detect the SNR required to generate 
this effect; we discuss this issue in §4.3 below. 

Recent evidence has suggested that characteristic ages can 
overestimate (Kaspi et al. 2001b) or underestimate (Gaensler & 
Erail 2000) a pulsar’s true age by up to an order of magnitude, 
depending on assumptions made about the initial period of and 
braking torque acting on the neutron star. The PWN / reverse 
shock interaction which we have proposed here is expected to 
take place ~ lO'* yrs after the supernova explosion (Reynolds 
& Chevalier 1984), in approximate agreement with the char¬ 
acteristic age r = P/2P = 21.4 kyr measured for PSR B1823- 
13. Thus within the reverse shock scenario, the available evi¬ 
dence suggests that the pulsar’s characteristic age and true age 
are broadly consistent. 

4.3. The Lack of an Associated Supernova Remnant 

The interpretation which we have adopted in §4.2 requires 
the presence of a surrounding SNR in order to generate the re¬ 
verse shock which consequently interacts with the PWN. How¬ 
ever, no such SNR has been detected at radio wavelengths 
(Braun, Goss, & Lyne 1989), nor do we see any surrounding 
shell of X-ray emission from such a source in the data presented 
here. 

The lack of detectable SNRs around pulsars with small char- 

* ’ http://asc.harvard.edu/toolkit/pimms.jsp 


acteristic ages has long been problematic (e.g. Reynolds & 
Aller 1985; Erail et al. 1995). It is thought that such SNRs 
must exist, as otherwise the PWN would expand freely into the 
ambient medium, and would produce a highly extended and un¬ 
observable “ghost nebula” (Blandford et al. 1973). In the spe¬ 
cific case of Vela-like pulsars such as PSR B1823-13, Braun 
et al. (1989) interpret the absence of a detectable radio SNR as 
resulting from the rapid expansion of the SNR blast wave into 
a low density cavity, followed by a collision with a dense sur¬ 
rounding shell evacuated by the progenitor wind. The SNR’s 
radio and X-ray emissions then fade rapidly as the SNR’s en¬ 
ergy is dissipated through radiative shocks. Simulations of this 
process clearly show that the SNR/shell interaction drives a re¬ 
verse shock back towards the center of the system (e.g. Tenorio- 
Tagle et al. 1990). Eurthermore, either an offset of the progen¬ 
itor from the center of its wind bubble or pressure variations 
within the turbulent SNR/bubble interaction region can cause 
different parts of the reverse shock to collide with the PWN at 
different times (Rozyczka et al. 1993; Dwarkadas 2001), com¬ 
pressing the PWN in a highly asymmetric fashion as is required 
here. 

We can consider whether such a SNR would be detected 
by scaling the emission seen from the Vela SNR. Lu & As- 
chenbach (2000) find that X-rays from the Vela SNR can be 
modeled as a two-component Raymond-Smith plasma, with 
temperatures for the two components of kTy « 0.15 keV and 
kT 2 ~ 1 keV. The volume-integrated emission measure for the 
cooler component is ~ 10 times that of the hotter compo¬ 
nent, and the total unabsorbed luminosity (0.1-2.5 keV) is 
2.2 X 10^^ erg s“* for a distance of 250 pc. Scaling to a distance 
of 4 kpc, the expected unabsorbed 0.1-2.5 keV flux densities 
for the two components are /i = 1.0 x 10“'° erg cm“^ s“' and 
/2 = 9.5 X 10“'^ erg cm“^ s“'. 

Using the W3PIMMS tool," we can determine that the ex¬ 
pected EPIC MOS count rate from such a source is 0.18- 
0.29 cts s“* for absorbing columns in the range (1.0- 1.4) x 
10^^ cm“^. Again scaling from the Vela SNR, we expect 
this emission to be distributed around a shell of diameter ~ 
20'. Assuming a shell thickness 20% of its radius, the ex¬ 
pected surface brightness for EPIC MOS would be (4-7) x 
10“^ cts s“' arcsec“^, several times below the observed back¬ 
ground level 10' off-axis of ~ 3 x 10“° cts s“' arcsec“^ (after 
correction for vignetting). We can make similar calculations 
for archival observations of this source with ROSAT and ASCA: 
for ROSAT PSPC we expect a SNR surface brightness (5 -9) x 
10“* cts s“' arcsec“^ against an observed background in these 
data of ~ 2 X 10“^ cts s“* arcsec“^, while for ASCA GIS we 
expect a SNR surface brightness (4-6) x 10“*' cts s“* arcsec“^ 
against the observed background of ~ 2 x 10“^ cts s“' arcsec“^. 
Thus we expect that any emission from a Vela-like SNR asso¬ 
ciated with PSR B1823-13 would not be detectable in any of 
these observations. This conclusion is strengthened by the fact 
that the characteristic age of PSR B1823-13 is approximately 
twice that of the Vela pulsar; if this indicates a similar differ¬ 
ence in these pulsars true ages, any associated SNR would most 
likely be even fainter than considered here. 

Eor the radio observations carried out in this direction by 
Braun et al. (1989), we estimate a 1-GHz surface bright¬ 
ness sensitivity for a SNR of this extent of S ^ l.Ox 
10“^* W m“^ Hz“* sr“' (3 a), where we have assumed a typ¬ 
ical radio photon index F = 1.5 and have accounted for pri- 
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mary beam attenuation. There are certainly many SNRs fainter 
than this, including relatively young SNRs such as G266.2-1.2 
(Duncan & Green 2000). Thus, the failure to detect a SNR asso¬ 
ciated with PSR B1823-13 is not inconsistent with the possibil¬ 
ity that a rapidly expanding, rapidly fading SNR has produced 
the pressure which has confined and distorted the PWN. 

5. CONCLUSIONS 

We have carried out a deep X-ray observation of the Vela-like 
pulsar B1823-13 with XMM-Newton. We find that all the de¬ 
tected X-ray emission results from the interaction between the 
relativistic pulsar wind and its environment: a bright elongated 
core immediately surrounding the pulsar most likely represents 
the wind termination shock, while a surrounding faint region 
with a softer spectrum corresponds to a pulsar wind nebula with 
a mean magnetic field of ~ 10 /xG. Emission from the pulsar it¬ 
self, either in the form of a central point source or as pulsations 
at the 101-ms pulsar period, is not detected. 

The nebula, which we designate G18.0-0.7, is distinctly one¬ 
sided: most of its emission is distributed in a region of extent 
5' restricted to the southern side of the pulsar. We interpret this 
morphology as resulting from the compression and distortion 
of the pulsar wind by the asymmetric reverse shock produced 
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by an unseen surrounding supernova remnant. Deeper obser¬ 
vations of PSR B1823-13 and other young pulsars at meter 
wavelengths, carried out with the Very Large Array and with 
LOFAR, are probably the most promising way to directly de¬ 
tect these associated SNRs. The reverse shock interaction be¬ 
tween a pulsar nebula and its surrounding supernova remnant 
is expected to occur 10-20 kyr after the supernova explosion, 
suggesting that many Vela-like pulsars might be in this phase of 
evolution. Future observations of other such pulsars and their 
nebulae with XMM-Newton can explore this possibility. 
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Table 1 

XMM-Newton OBSERVATIONS OF PSR B1823-13. 


Date of Observation 

Length of Observation (ks) 

Effective Exposure (ks) 
MOSl MOS2 pn 

2001 October 16 

30.5 

17.8 

18.3 

10.2 

2001 October 18 

29.8 


23.2 

12.6 

TOTAL 

60.3 


41.5 

22.8 


Note that the effective exposures have been corrected for dead-time effects. 


Table 2 

Radio timing ephemeris for PSR B1823-13. 


RA (J2000) 

Dec. (J2000) 

Erequency (Hz) 

Erequency derivative (Hz s“') 
Erequency second derivative (Hz s“^) 
Epoch of Ephemeris (MJD) 

Range of Validity (MJD) 


18''26“13?175(2) 
-13°34'46"7(2) 
9.85564535524(1) 
-7.2930848(2) x IQ-'^ 
1.51530(6) X 10-22 
50750.000000 
49156 to 52359 


These data have been supplied by M. Kramer from timing observations with Jodrell Bank Observatory. Numbers in parentheses indicate the uncertainty in the last 
significant figure of each value. 


Tables 

Spectral eits to X-ray emission erom G18.0-0.7. 


Region 


Total Counts 


Model 

N„ 

T/kT/EQ 

fx 

xi/v 


MOSl 

MOS2 

pn 


(10^^ cm“^) 

—/(keV)/(keV) 

(I0“” ergcm^^ 


Core^ 

455 ± 23 

458 ± 22 

8I5±33 

PL 

L0±0.2 

1 A+O.* 

L6_q2 

3 5"’’' .'’ 

^..^-0.8 

42/38= 1.1 





PL 

1.0 (fixed) 

I.6±0.I 

2.9 ±0.3 

42/39= 1.1 





PL 

1.4 (fixed) 

I.8±0.I 

3.0 ±0.3 

49/39= 1.3 





BB 

< 0.2 

l.I9±0.08 

2.1 ±0.5 

60/38 = 1.6 

Diffuse 

1843 ± 135 

2360 ± 138 


PL 

1 4+''-^ 
^•^-0.2 

'J-I-U.4 

2^-^-0.3 


130/133= 1.0 

Component 




PL 

1.0 (fixed) 

I.9±0.1 

13±2 

135/134= 1.0 





PL 

1.4 (fixed) 

'y 'i-l-0.1 

2^-^-0.2 

17 ±3 

130/134= 1.0 





Brem 

^•^-0.2 

3.9!?;? 

12+^ 

129/133 = 1.0 





SRCUT 

1.4 ±0.2 

~ 1 

I5±I 

130/133 = 1.0 





RS 

1 

_-0.2_ 

3 4+1-3 

133 

133/133 = 1.0 





MEKAL 


133 

136/133 = 1.0 





NEI 

__ 

383 

126/132= 1.0 


Uncertainties are all at 90% confidence. All models assume interstellar absorption using the cross-sections of Balucihska-Church & McCammon (1992), assuming 
solar abundances. Models used: “PL” indicates a power law of the form fe oc where T is the photon index; “BB” indicates a blackbody emitting at a temperature 
T ; “Brem” indicates thermal bremsstrahlung; “SRCUT” indicates a synchrotron spectrum from a power-law distribution of electrons in a homogeneous magnetic 
field, with an exponential cut-off at a characteristic energy Eq (Reynolds & Keohane 1999); “RS” indicates a Raymond-Smith spectrum of temperature T (Raymond & 
Smith 1977); “MEKAL” indicates a hot plasma model at temperature T ; “NEI” indicates a non-equilibrium ionization plasma at temperature T with solar abundances. 

"Flux densities are for the energy range 0.5-10 keV, and have been corrected for interstellar absorption. 

^Counts and flux densities listed here for the core correspond only to 75% of the total for this region (see text for details). 
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18’^26™30® 20® 10® 26™0® 50® 


RA (J2000) 

Fig. 1.— EPIC MOS image of the field surrounding PSR B1823-13, produced by combining all MOSl and MOS2 data listed in Table 1 in the energy range 
0.5-10 keV. The image has been corrected for off-axis vignetting, and then convolved with a gaussian of FWHM 20''. The image is displayed using a linear transfer 
function ranging between 0% and 30% of the peak intensity in the field; the white contours are plotted at levels of 13%, 15%, 20%, 35%, 75% and 90% of the peak. 
The white cross marks the position of the pulsar as listed in Table 2, while the circle shows the outer edge of the extraction region used for the spectrum of the 
diffuse component of the source. Linear structures in the image correspond to gaps between individual CCDs. 
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RA (J2000) 

Fig. 2.— As in Figure 1, but showing emission from the central 90" x 90" region surrounding the pulsar, and smoothed with a gaussian of FWHM The 
image shows a linear transfer function ranging between 0% and 100% of the peak intensity. The white cross marks the position of the pulsar as listed in Table 2, 
while the circle shows the extraction region for the spectrum of the core. 
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-20 0 20 
Offset (aresee) 

Fig. 3.— Intensity profiles through the core region of G18.0-0.7. The histograms plot the summed flux of EPIC MOS data in successive binned regions along 
each of the RA and Dec axes. The size of each region is two MOS pixels (2F2) along the axis of offset, and three MOS pixels (3F3) along the axis for which the data 
are summed. Positive offsets along the RA axis (upper panel) correspond to more westerly points, while positive offsets along the Declination axis (lower panel) 
correspond to more northerly data; in both cases the pulsar position is at zero offset. The smooth curve on each plot is the profile of the EPIC MOS point spread 
function at an energy of 1.5 keV, after similar summing as applied to the image. 
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-200 0 200 
Offset (aresee) 

Fig. 4.— Intensity profiles through G18.0-0.7 along the north-south direction. The solid histogram plots the summed flux of EPIC MOS data in successive binned 
regions along the Declination axis. The size of each region is 16 pixels (17^6) along the offset axis, and 208 pixels (3f8) along the axis for which the data are 
summed. Increasing offsets correspond to more northerly data; the pulsar position is at zero offset. The broken line corresponds to an identical analysis for blank 
field observations, after appropriate normalization (see text for details). 
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